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Sixteen points are located on a 100 x 100 grid, 
as shown on the cover. The points are to move as follows: & 
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Each point moves 1/7 of the way toward the next higher 
numbered point. Thus, for example, point No. 5 should 
move to a new position given by: 


bed 
i] 


1 
= 58 + (75 - 58) 


y = 22 oe - 22) 


which puts it at (60.42857, 26.14286). Point No. 16 
moves in the direction of point No. 1. The set of points 


is to be moved in order, from No. 1 through No. 16, with 
the process then repeating, starting again with point No. l. 


PROBLEM 159 


It seems intuitively clear that the process will 
converge; that is, all the points will eventually meet. 


If so, where? 
Will it be close to the center of the initial set of 


16 points; that is, at the average of the initial x and y 
values, which is 55.3125 and 48.6875? 
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An Approach to Floyd's Problem 
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Last month I discussed some of the virtues of Floyd's 
problem, which is to divide the numbers 


0, Vey oo VO 


into two parts whose sums are as equal as you can make 

them after ten seconds of computer time. Here is the way 
I decided to approach the subject; perhaps some reader will 
have a much better idea. 


Since (V1 +VY2+ ... +¥V50) is equal to 


119.51790 03017 60392 24702 02231... , 


we seek the subset ye Rib, VER ... » ¥50} whose sum least 
exceeds this value. 


In the first place it is helpful to estimate the kind 
of results that might be expected. Most subsets of the 
given set have about 25 elements, and in fact the number 
of subsets with exactly 25 elements is 50!/25!25!; this is 
about 


200 
257 


according to Stirling's approximation. [The exact number 
4s 12641 06064 37752, compared to 


3 


299 = 1 12589 99068 42624] 


All of these subsets have a sum which lies between 


VelaVeg ... 25.=985.6 anay V26 +74 22 4VEo anisole 


and they tend to cluster about the average value 119.5. 
So we have more than 1014 numbers packed into an interval 


of length less than 70; therefore 70/1014 appears to be 
a conservative estimate for the amount by which the best 
partitioning exceeds 


ae PY ee hs +50). > 
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But not all of these subsets will give different 
sums, In the first place, the seven numbers 


Wile ee VAG = 152°. ..,7 


will yield only integer values (and it is easy to verify 
that each integer from O to 28 can be represented). Our 
problem therefore reduces to finding a subset of the 43 


numbers 
(V2. V3. V5> Ve, ... » Va8, V50} 


whose sum has a fraction part least exceeding 
51790 03017 60392 24702; 


unless the integer part of the sum is unusually large or 
small, we will be able to adjust it by addins: a suitable 
subset of {1,2,...,7}, bringing the integer part up to 119. 


We now have 243 possibilities to try; but they aren't 
all distinct, either. maine values 


\/2, V8, Vi8, V32, V50} = {Ve. aVe, 3V2, 4Ve sve} 
lead to subsets whose sum is restricted to sixteen values 


Oe 2 Yen -. 152, 


so only 16 of the 32 subsets are different. Similarly, 
there are only eleven essentially different subsets of 


V3, V12, V27, V48 


to try. By means of these observations our program can 
do in 10 seconds what it would take (10:32-16)/(16°-11) * 29 
seconds to do otherwise. 


But there are still more than 241s 10914 possibilities 
remaining, and this is far too large; if it takes 100us 5 
for me to test one subset, I will have time to test only 10 
subsets. 


One way to gain speed is to divide 


VE, V5, ... » Vi8, V0} 


into two parts A and B, and to form a table containing 
all sums of the subsets of A. Then for each subset of B 
we can look in the table to see if there is an entry with 
the appropriate leading bits. 


That is what I eventually did. Let 
a =/-V3, Ve, Vib, V32, V50, V3, Vie, Ver, 
Via, V5, V6, V7, Vio, Vin, V3, Vis, Vis}, 


so that the subsets of A take on (11/16) -226 different 
values. I stored the 32 leading bits of these fraction 
parts into a hash table of size alo, and in another table of 
size 216 I stored a bit pattern identifying the corresponding 
subset. (The hash code was the low order 16 bits of each 


32 bit key.) Let B ={V17, Vi9, V20, ... , V47} be the 


remaining set of 26 elements; in the machine I used only , 
the 32 leading bits of the fraction parts of each number Vk. 


Since my hash table contained (11/16) +226 more or less 
random 32-bit numbers, I expected to get a match once out 


of every 232/( 11/16) 216 = (16/11) -27° times I looked up 
another 32-bit value. Therefore I wanted to have a program 
that generated 2l7 or so subsets S of B, looking up the 32 
bits corresponding to (.51790 ... - =S)mod 1 in the hash 
table. I figured that would probably give me two or three 
matches, and I could choose the best corresponding partition. 
But I would have only about 75 to 80 us to spend per subset, 
so I needed a fast way to probe the hash table and to 
generate the subset sums. 


The solution was to use ordered hash tables with linear 
probing; this is a variant of hashing which Ole Amble and I 
had discovered a few years ago {2 : Such hash tables are 
especially well suited to cases when searches are unsuccess- 
ful, requiring only 2.1 probes per search in this case. 

For the subset generation, I used Gray code fell since this 
meant that only one subset element changed state each time 
and the subset sum was therefore easy to update. It 

took about 1.5 seconds to build the hash table. I started 
the subset generation of B in a random part of its cycle, 
since I knew that I would only be able to look at a small 
fraction of its subsets. 


The best of the three results I got before 10 seconds 
expired was actually typed out after only 6 seconds, namely 


V2 +V3 +Vi +V¥5 +V6 +V12 +-V17 +-V19 +V22 +V27 + 
V28 + -V29 +V33 + V34 + V35 +V37 +V38 +V49 + V50 


== 119.51790 03021 65123 39726 54768 


This sum differs from its complementary sum by approximately 
8-10729, (The optimum partition might well be a thousand 


times better, as mentioned above; my program could have 
found it if 10 thousand seconds were available instead of 
ten!) 
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I used a fairly large and fast computer, the DEC KL-10 
at Stanford's Artificial Intelligence Laboratory. ate aC 
had used a slower machine, I would have cut down the number 
of bits in each hash table entry from 32 to something less, 
one bic per factor of two in speed. If I had used a 
smaller machine, I would have had to make the hash table 
smaller. 


If I really wanted to get the optimum partition, I 
think it could be found in about 30 minutes. The best 
approach I can think of would be to divide the 43 irrational 


numbers into two subsets having respectively 1/16 and get 


distinct sums. I would sort each of these files of sums 
by their fraction parts; since each file contains about two 
million entries, this would be the time-consuming part of 
the operation. (Instead of using a general-purpose sort 
routine I would write a special one, since it is easy to 
go from the sorted file for set A to the twice-as-long 
sorted file for Auja} by essentially merging the first file 
with itself. The total time to build the sorted set for A 
when A has n elements is therefore of order 2", while a 
general-purpose sort routine would take order n-2" steps.) 
Finally, given two large sorted files 


xy iene sect Xx and Vy & oho SH 


hal n 


and a number z, there is a nice algorithm which computes 
min{x, + yy | xy + Y,2 24. 
The reader will enjoy discovering this algorithm for 
himself; curiously it is essentially the same as Hamming's 


"p/q" algorithm in PC41-4, under a logarithmic transform- 
ation, with x, = log p and Oe log(ax). 


I tested my hashing program by trying it first on 
Wa Veh «8 50} 
and then on {Vi, V2, Eckee 35 Vito}. 


I gave the optimum partition for the former case in the 
previous article (see issue No. 46); and I think I found 
the optimum partition also in the latter case: 


Ve +5 +V7 +V13 +V14 + VI7 +V18 + Vig +V¥22 +V23 + 
V24 +V26 + V27 +V29 +30 +V32 +V34 +V38 +V39 


= 85.807894002346 , 


Vi +-V¥3 +-V4 +V6 +V8 404... +ViI0 = 85.807894000154 


However, I am not absolutely sure that this is the best, 
because there is a slight chance that an unusual rounding 
error might have occurred somewhere in the calculations. 
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I still think that the 10-second-limited problem is 
more interesting; according to Hamming's famous aphorism, 
the purpose of computing is insight, and I believe this 
problem brings some valuable insights into view. 


Postscript: R. L. Graham recently told me about the 
following sets of nine square roots whose sums agree to 
54 decimal places: 


¥100000001 +7°¥100000025 +7¥100000031 +V¥100000084 + 
100000087 +7°7100000134 +¥100000158 +V100000182 + 
V¥100000198 
= 9000 .00449 99983 56751 33987 35593 13014 12703 30519 
82221 56985 62024 86633 


100000002 +°V¥100000018 +Vi00000042 +°¥100000066 + 
7100000113 +¥100000116 +¥100000169 +V¥100000175 + 
¥100000199 
= 9000.00449 99983 56751 33987 35593 13014 12703 30519 
82221 56985 62028 18260 
To find the secret of his construction, see Hardy and 


Wright's Theory of Numbers, 4th edition, p. 338 (notes on 
Section 21.10). 
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2-3-5 again 


In issue No. 10 there first appeared this problem: 


Show the logic of generating, in 
order, all numbers having only the 
factors 2, 3, and/or 5. The first 
30 such numbers are: 


PS a 50300, Gyr 10h, LA ee. Te. 
18,20, 24, 25, 27, 30) 32) 36. 40% 
45, 48, 50, 54, 60, 64, 72, 75, 80, 81. 


The problem was attributed to Richard Hamming. 


In issue No. 16, three different solutions to the 
2-3-5 problem were shown, including one by Hamming. The & 
following comes from Norman Sanders, Trondheim, Norway. 


I had the problem of generating thousands of numbers 
of the form 


Sone peigh tae: 28.3? .5°.74 


in connection with the numerical validation of an 
approximation in prime-number theory. But that was in 

the mid-50's on EDSAC-1, which had a 1K store, no tapes, 

no divide, 1.5 millisecond cycle time, and zero-generation 
reliability. A divide subroutine was a time-consuming 
luxury that ruled out the first two solutions shown in 
PC16. The severe storage limitations ruled out long lists, 
and the low availability of the computer implied a simple 
program. 


As an example of the general method, the solution for 
go Re 


was to generate three lists, A, B, and C (of which C 

consisted of a single entry only) and to transfer elements & 
from list to list, freeing space at each transfer and 

each print instruction. In this way, I was able to 


exceed 210° without reaching the limitation of the store. e 
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Snapshot at the time of printing N = e5: 


PRINT LIST A-LIST B-LIST C-LIST 
2 0.0 010 eiptoreal 
3 Mees @ a2 
4 00 Onane 00 3 
5 oe pt oom 
6 sigalg Duo 2 
8 300 30 
9 2u0 Ge 1 
ol 
1 


et fj 
Oo F OO OO fF O 


ot, 


10) 

10 1 012 

12 2 o40 

15 Ol 00 3 

16 4 0 © 5) sl 

18 2 € 
20 0) 

au 1 

25 0 


=- OO WwW YD WN FY OO 
Pewee Or © mM Ww 
© [h(i oe an? 


values are 22.30.5¢ in increasing order of value. The 


The A-list consists of triples (a, b, c), whose 
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B-list consists of triples (0, B, 8 ), whose values are 


Q 
close in order. The C-list contains the single triple 


(0,0,°), of value 5, Initially the A-list consists of 


the single entry (1, 0, 0); the B-list of the single entry 


(0, 1, 


0), and C of the entry (0, 0, 1). The logic of 


the algorithm is shown in the flowchart. 


Te entries in the lists above the dotted lines have been 


removed. N = 25(0,0,2) has just been printed, and the 


next step is to compare (1,0,2) (= 50) with (0,1,2) (= 75) 


at the top of the B-list. It is less, so (1,0,2) will 


be appended below (4,1,0) in the A-list. 


Intrusion, 


A few pertinent remarks: 


Don't throw anything away. I'm far too 
experienced today to reproduce such a 
simple method. 


Document your stuff so that you can 
understand it twenty years later. You 
never know what the Hamming birds are 
going to come up with. 


When we had hardly any computers, the 
algorithms we used had to treat the 
idlosyncracies of the hardware very 
carefully. Now that we have thousands 
of computers to choose from, most 
programming seems to be machine- 
independent. Computing has become 
Aristotelian. 


We managed to get a lot of useful work out 
of the primitive computers. If we hadn't 
they would have stayed primitive, or died 
out all together. aa 


McGraw-Hill, 1973. His pixie-ish lack of respect for authority should 


not detract from his computing wisdom. 
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In Martin Gardner's "Mathematical Games" column in the 
November 1976 Scientific American, he shows that the following 
trick works: 


Ask someone to pick a number, write it down, and show it to 
you. You add one to his number and divide by some other number 
that you have chosen (and secretly written down) and give him the 
result. Your friend now has two numbers and is to follow this 
plan: Given two numbers, add one to the second number, divide by 
the first number, and record the result as the third number. 
Then, using the last two of these, repeat the procedure two more 
times; the final result will be your hidden number. 


The arithmetic for the choice of 17 by the victim 
and 1.2345679 as your hidden number is shown here: 
Q) Your hidden number: 1.2345679 
@) His first number: 17 


@) You give him ss Dane 

@) re = .9164705882 
© Ser ei = .1314451707 
© SE = 1.2345679 


It is easy to show that this must work. The 
successive stages for any two starting numbers will be: 


1 
@) y : (a) ytxrel 
() — (6) oe 


proBLem 160 


Thus, the numbers will repeat on a cycle of 5. 


Problem: To find a simple similar algorithm that lead 
will repeat on a cycle of 6. 


Problem Solution 
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@ Wendy's Problem (Number 148, PC45-6) called for 
counting all possible triangles that can be formed on 
a 10 x 10 grid of points and tabulate the results by 


areas. 


Associate Editor David Babcock produced the 
tabulation shown here. The column labelled N is twice 
the area and column K shows how many such triangles 
there are. The total of the K values is 161,700. 


It is of interest to note that there are 4448 
combinations of 3 points that do not form a triangle, 
and the curious distribution of the other 80 possibilities. 
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F-N SEQUENCES 


In the sequences shown in Table A, F-1 is the 
familiar Fibonacci sequence, in which 


es = Fy + Fo 


It is well known that the ratio of one term to the 
preceding term 


FL /Fyy 


converges to ine = 1.618033988... 


The sequence F-2, starting with three ones, is formed 
by the relation 


FE = Fyo + Fy_3 


and it appears to converge to the ratio 1.3247... 

We say “appears to converge" because the ratio of one 

term to the preceding term has a cycle of two values. For 
example, successive ratios for the first few terms are: 


Hepiay Ss. 1), 3333 
21/16 = 1.3125 
28/21 = 1.3333 
Br /ea ss) 21 3214 
49/37 = 1.3243 
65/49 = 1.3265 
66/65 = 1.3231 


114/86 = 1.3256 

151/114 = 1.3246 

200/151 = 1.3245 e@ 
so that the value to which the ratio converges is the 


average of two ratios. Similarly, for F-3, the convergent » 
value will be found by averaging three successive ratios. 
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expresses the relationship shown in the curve? Is the 
curve asymptotic to y = 1? 


N-senies 47 


Subsequent members of the F-N family of sequences 
are defined similarly: 


F-3: 
F-4; 
F-5: F 
F-6: F 


Fy-3 a he ih converging to 1.2207.. 


tt 


Rael Byes converging to 1.1672... 


eetectskn G converging to 1.1347... 


Fy-6 + Fn-7 converging to 1.1128... 
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For the sequences F-1 through F-8, the convergent 
values plot as shown. What is the curve that best 


i 


i 
LH) F 


16720978579 35717464414219 3994492006401598030984 29948 
3 .85014760171005858682095066977 217 3708896050502020224 
6.8556546004010441249 3587 144908484 896046064 3461001326 
3 .6088260801 38694689252517295958892614905551690162338 
1.4696 36013359929253623112224 37 1627488588820827 562794 
1..0392522623109629533847966292221986501217 36878483157 


2581312886190067 39623 . 28580021527 33804 3163708299 30440 
608106139724 368227 5040978221987 


232297 22223604.1657886 385 .7994890281440584469844961463 
1767505405499833821420635585 


1.5495229407708354.15559197706584 meena i 


Cunningham's Process 


In the May, 1962 American Mathematical MONTHLY, 
G. S. Cunningham showed that any positive fraction can be 
represented as a sum cf distinct unit fractions. 


a 1 1 1 lt 

ee ee ee 

b b b b b 
a terms 


Retain one of the duplicating terms. Replace each of the 
others by the identity 


1 al 1 


—_ Bee + ——= 
n nt+1 (n+1)n 
and repeat this process until all terms are distinct. 


For example: 


aie Es eee 
o_ = — + — 
3 3 3 
fee 
3. 2 Wie 
321,1,2 
4 4 4 4 
ot Poe ee 
TONS en Tees aaa) 
a eka eel 1 1 1 
piwiaee 0% G0 40:) |. aT alan 
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1 

Ps ee BO wom SO GSO 

eee ee ee te 

GemiG cor .7 te 31° 930° 7 * Ge * 31 * S30 

eee ee  . atk SS 
Mimecoeeses 7 42 31.930 8 56° 43” Uaede 
mew 1 1 2 


+ ee 
oe SYe°S Isha 931°930 


Se 4 et ee 
fe «31° 7930 6. 560 22 1806 


ik Ba) es Ht i 
32 gg2. 931 865830 


WN 
In each example, the fraction used was < ‘0 
= 
Problem: extend this table: a 
© 
Number of Largest o. 
N final terms denominator 
3 Bs] 12 
4 vi 420 


5 15 865830 
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